Segregated tunneling-percolation model for transport nonuniversality 
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We propose a theory of the origin of transport nonuniversality in disordered insulating-conducting 
compounds based on the interplay between microstructure and tunneling processes between metallic 
grains dispersed in the insulating host. We show that if the metallic phase is arranged in quasi-one 
dimensional chains of conducting grains, then the distribution function of the chain conductivities 
g has a power-law divergence for y — > leading to nonuniversal values of the transport critical 
exponent t. We evaluate the critical exponent t by Monte Carlo calculations on a cubic lattice and 
show that our model can describe universal as well nonuniversal behavior of transport depending 
on the value of few microstructural parameters. Such segregated tunneUng-percolation model can 
describe the microstructure of a quite vast class of materials known as thick-film resistors which 
display universal or nonuniversal values of t depending on the composition. 
PACS numbers: 72.60.-hg, 64.60. Fr, 72.80.Tm 



I. INTRODUCTION 

When the conductivity cr of an insulating-conducting 
compound is measured as a function of the volume con- 
centration p of the conducting phase, one finds that by 
reducing p eventually the system undergoes a conductor- 
to-insulator transition at a particular critical value pc 
of the volume concentration. In the critical region < 
p — Pc"^ 1 the coirductivity follows a power law behavior 
of the form: 



p{9)^pK9) + {^~p)5{9), 



(2) 



where 5{g) is the Dirac delta-function and h{g) is the 
distribution function of the finite bond conductivities. 
For well behaved h[g), transport is universal and follows 
Eq.(l) with t = to ~ 2.0 for three dimensional lattices. 
Instead, if h{g) has a power law divergence for small g of 
the form: 
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where ctq is a prefactor which depends on the particular 
system considered aird i is a positive number typically 
larger than the unity. 

Percolation theory explains the power law form of 
Eq.(l) as being due to the lack of any cut-off length 
scale apart the linear size of the sample and predicts 
that the exponent t is universal and depends only upon 
the dimensionality of the system.^ This prediction is con- 
firmed by various granular metals compounds and model 
systems which have been found to follow Eq.(l) with 
t ~ 2.0,^"^ that is the value obtained by numerical cal- 
culations on three-dimensional random resistor network 
(RRN) models.^ 

In addition to systems showing universality, a large 
number of disordered compounds displaying values of t 
larger than t ~ 2.0 have been repeatedly reported, ''"^'^ so 
that in the present situation it appears that t can assume 
any value between t ~ 2.0 up to about t ~ 6.0 — 7.0. 

Within percolation theory on a RRN, Kogut and 
Straley showed that universality breakdown of trans- 
port may arise from anomalous distributions of elemen- 
tal conductivities.^"* By assigning to each neighbouring 
couple of sites on a regular lattice a bond with finite 
conductivity g with probability p and zero conductivity 
with probability 1 — p, the resulting bond conductivity 
distribution function becomes: 



and a is larger than a critical value Uc, Kogut and Straley 
showed that transport is no longer universal and the con- 
ductivity exponent becomes dependent on a}^ Renor- 
malization group analysis predicts in fact that 
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where D is the dimensionality of the lattice and is the 
correlation-length exponent {i^ — 4/3 for D = 2 and ly ~ 
0.88 for D = 3).i5"*'^ For £» = 3 and by using - 2.0 
and ~ 0.88 the critical value of the exponent is ac — 
0.107. 

Microscopic models which may justify Eq.(3) are the 
random void (RV) model proposed by Halpcrin, Seng and 
Fen,*^ and the tunneling-percolating model of Balberg.*^ 
The RV model describes a system of insulating spheres 
(or disks in two dimensions) embedded randomly in a 
continuous conducting material. In this situation, trans- 
port is dominated by the conductivity of the narrow necks 
bounded by three interpenetrating insulating spheres. 
Such necks have a wide distribution in widths resulting 
in a wide distribution of conductivities. The original for- 
mulation of the RV model predicted t = to + 0.5 for the 
conductivity exponent of the whole sample. A recent gen- 
eralization of the RV model by Balberg has shown that 
t can assume even higher values aird that in principle is 
not bounded above. 
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In the tunncling-pcrcolating model of Ref. [ 11], trans- 
port is assumed to be dominated by quantum tunnel- 
ing between neighbouring conducting particles dispersed 
in an insulating medium. If the distribution function 
P(r) of the distance r between two neighbouring parti- 
cles decays with r much slower than the tunneling decay 
exp(— 2r/^), where £, is the localization length, then the 
tunneling conductivity distribution function h{g) can be 
shown to behave as Eq.(3) with a ~ 1 — ^/2a, so that 
the transport exponent t becomes dependent of the mean 
tunneling distance o.^^ Interactions between the conduct- 
ing and insulating phases as well as properties of the mi- 
crostructure are argued to concur to the r-depcndence of 
P{r). Due to the complexity of the problem, explicit cal- 
culations of the interparticle distance distribution func- 
tion are missing and one must relay on phenomenological 
forms of P{r). 

In this paper we provide a microscopic derivation 
of -P(r) which has been inspired by the peculiar mi- 
crostructure observed in a particular class of insulating- 
conducting compounds: the so-called thick-film resis- 
tors (TFRs). These compounds are based on Ru02 (or 
Bi2Ru207, Pb2Ru206, and Ir02) grains mixed and fired 
with glass powders.^" Typically, TFRs are often in a seg- 
regated structure regime in which large regions of glass 
constraint the much smaller conducting grains to be seg- 
regated in between the interstices of neighbouring glass 
grains. Micrographs reveal that the conducting grains 
are arranged in a network of filaments spanning the en- 
tire sample. ^'^^'^^ By taking into account the quasi-one 
dimensional structure of such filaments and by neglecting 
interactions with the insulating phase, we show that the 
resulting P(r) can decay much slower than the tunneling 
decay leading to nonuniversal behavior of transport. 

This paper is organized as follows. In the next section 
we construct a RRN model which captures the essential 
structure of the filamentary network of TFRs and cal- 
culate the resulting P(r) and the distribution function 
of the conductivity of filaments. In Sec. Ill we perform 
Monte Carlo calculations and calculate the conductivity 
exponent t for a variety of situations. The last section is 
devoted to discussions and conclusions. 



II. THE MODEL 

Before describing our model in details, we find it useful 

to first discuss in general the interplay between the spa- 
tial distribution of the conducting phase within the insu- 
lating matrix and transport properties. Let us consider a 
generic insiilating-conducting compound where the con- 
ducting grains are embedded in an insulating host. In 
this situation, electron transfer is governed by electron 
tunneling from grain to grain. The grain charging energy 
and the Coulomb interaction between charged grains af- 
fects the overall transport properties especially regarding 
their behavior in temperature. Here we focus on systems 



where the temperature is high enough to possibly neglect 
charging and Coulomb effects so that the main electron 
transfer is dominated solely by tunneling leading to in- 
tergrain conductivity of the form: 

a{r) = aoe-2('-*)/«, (5) 

where (Tq is a constant which can be set equal to the unity 
without loss of generality, ^ oc is the tunneling 

factor (or localization length) and V is the intergrain 
barrier potential. In Eq.(5) we have approximated the 
conducting grains by spheres of diameter $ and r is the 
distance between the centers of two spheres which we 
treat as impenetrable (r > <&). 

Due to the exponential decay of Eq.(5), contributions 
to cr(r) from far away spheres can be neglected, -^^'^^ 
so that from now on r denotes the distance of two 
nearest-neighbouring spheres. Hence, the ensemble de- 
pendence of a{r) upon r is completely defined by the dis- 
tribution function P{r) of the distance between nearest- 
neighbouring spheres. In fact, once P{r) is known, the 
conductivity distribution function h{g) can be obtained 
as follows: 

h{g)= j drP{r)5\g - u{r)]. (6) 

In this preliminary discussion, we are interested in study- 
ing how the form of P(r) affects h{g) via Eq.(6) and 
which arc the requisites of P{r) which (;vcntually could 
generate a power-law distribution function as that of 
Eq.(3). As already pointed out, P(r) depends on the 
microstructure of the composite and on eventual interac- 
tions between the insulating and the conducting phases. 
In principle, therefore, the form of P{r) depends on the 
particular composite considered. However, if we imag- 
ine that interactions can be neglected, then it is natural 
to assume that the conducting spheres arc Poisson dis- 
tributed within the insulating phase. Then if D is the 
dimensionality of the system, by following Refs. [ 24, 25] 
the nearest-neighbour distance distribution function is 
approximatively of the form: 

P{r) , (7) 

where au is a constant depending on the mean distance 
between neighbouring spheres. Equation (7) is an asymp- 
totic approximation of the true P(r) and is valid only in 
the r/ao 3> 1 limit. This is however the limiting region 
of interest to us since it governs, via Eci.(5), the 
regime. It is also worth to point out that Eq.(7) holds 
true for penetrable as well as impenetrable (hard-core) 
spheres, the only difference being in the; explicit expres- 
sion for ao which is however of not importance at the 
moment. 

By inserting Eqs.(5,7) into Eq.(6), the resulting con- 
ductivity distribution function becomes: 
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which after some manipulations reduces to: 
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(9) 



For D = 2 and D = 3, the g ^ limit of the above 
expression goes to zero irrespectively of the value of 
^/2aD. In this case therefore no power-law divergence 
of h{g) is encountered and, as discussed in the introduc- 
tion, transport is governed by the universal critical ex- 
ponent t = to c:^ 2.0. Instead, when D = 1, equation (9) 
becomes: 
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which is exactly of the form of Eq.(3) if we identify a 
with 1 — £,/2ai. We have arrived therefore at the result 
that if the spheres are Poisson distributed along a one- 
dimensional line, the resulting conductivity distribution 
function has a power-law behavior for small g and, con- 
sequently, transport is nonuniversal for sufficiently large 
values of 1 — £,/2a\. 

The difference between the -D = 2, 3 and the D = 1 
cases stems from the decay of Eq.(7) which for _D = 2, 3 is 
much too fast with respect to the simple exponential de- 
cay of Eq. (5) . In fact from Eq. (6) it is simple to show that 
as long as Imir^oo Pir) / <y{r) = then limg^o ^(5) = 
irrespectively of the detailed structure of P(r). Hence to 
construct a RRN model having h{g) of the form of Eq.(3) 
we must consider forms of P(r) whose decay for r — > 00 is 
sufficiently slow. The result of Eq.(lO) suggests that for 
this scope one-dimensionality is an important ingredient, 
at least as long as interactions between conducting and 
insulating phases can be neglected. 

Among the various insulating-conducting compounds, 
thick film resistors are systems whose microstructure 
can be appropriately described in terms of quasi one- 
dimensional units. Let us consider the highly non- 
homogeneous microstructure typical of TFRs. These sys- 
tems are constitutcid by a mixture of large glassy parti- 
cles (with size L of order 1-3 jim) and small conduct- 
ing grains of size $ typically varying between ~ 10 nm 
up to ^ 200 nm. Due to the high values of L/<^, the 
small metallic grains tend to occupy the narrow regions 
between the much larger insulating zones leading to a 
filamentary distribution of the conducting phase. ^'^-'^'^^ 
A classical model to describe such a segregation effect 
was proposed already in the 1970's by Pike.^ This model 
replaces the glassy particles by insulating cubes of size 
L ^ <I> whose edges can be occupied by chains formed 
by adjacent metallic spheres of diameter Let us as- 
sume that an edge has probability p of being occupied 




(b) 



FIG. 1. (a): pictorial representation of the segregated tun- 
neling-percolation model. The cubes represent insulating 
grains while the spheres are conducting particles. The spheres 
are arranged to occupy the edges of the insulating cubes with 
probability p. The total ensemble of occupied an unoccu- 
pied edges forms a cubic lattice spanning the entire sample, 
(b): equivalence between an edge occupied by n -|- 1 spheres 
and a conducting element. The set of inter-spheres tunnel- 
ing conductivities is equivalent to a conductor with n resis- 
tive elements with conductivities cr, in series. The fluctuation 
in distance between two neighbouring spheres leads to fluc- 
tuating tunneling conductivities. With this equivalence, the 
model depicted in (a) can be considered as a bond-percolation 
model whore a fraction p of bonds has variable conductivities 
and a fraction 1 — p is insulating. 



by a chain of n -I- 1 spheres and probability 1 — p of be- 
ing empty. As depicted in Fig.l, the set of occupied and 
empty edges form a cubic lattice spanning the entire sam- 
ple. 

To define the RRN relevant for this model we pro- 
ceed as follows. The conductivity of a single occupied 
channel is governed by the conductivities of the metal- 
lic spheres and those between pairs of two-neighbouring 
spheres. The conductivity of the metallic spheres adds 
only a negligible contribution to g which is then given 
by n conductivities o-j of pairs of nearest-neighbouring 
spheres in series: 



g = 



n ^ 



(11) 



We assume that the inter-sphere conductivities crj are due 
to tunneling processes between two adjacent spheres so 
that their sphere-to-sphere distance r dependence is that 
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of Eq.(5). For TFRs, the tunneling hypothesis is well 
sustained by their high values of piezoresistance (i.e., 
the strain sensitivity of transport), and the low tem- 
perature dependence of transport indicating some kind 
of assisted hopping. As done in the introductory part of 
this section, we neglect interactions between the insulat- 
ing and conducting phases, and assume that the sphere 
centers are Poisson distributed along the cube edge. In 
doing so, wc implicitly assume that finite size effects of 
the channels can be neglected and that periodic bound- 
ary conditions are applied. In this way the last sphere on 
one end of the channel is identified with the first one on 
the opposite end, so that we have n individual spheres 
and n inter-sphere tunneling junctions. In this situa- 
tion, the distances r change according to the distribution 
function Pn{r) of the nearest-neighbour distances r of n 
impenetrable spheres arranged randomly in a quasi one 
dimensional channel. By following Ref.^^, Pn{r) can be 
calculated exactly and it is given by: 

P„(r) = — ?_e-('-*)/('^"-*)e(r - $), (12) 
a„ — $ 



where 9 is the step function and 



(13) 



is the mean inter-sphere (center-to-center) distance. In 

the above expression n^/L cannot be larger than the 
unity since no more than L/$ spheres can be accommo- 
dated inside a channel. Note that the asymptotic ex- 
pression Eq.(7) for D = 1 coincides with Eq.(12) if ai is 
identified with a„ — $. Hence the distribution function 
/(cr) of the intcr-sphcre conductivities should be of the 
same form of Eq.(lO). In fact: 

poo 

f{a)= drP„(r)%-a(r)] = (l-a„)a-"", (14) 
Jo 



where 



e/2 



(15) 



Having obtained an explicit expression for the distribu- 
tion function /(cr) of the inter-spheres conductivities, we 

can now calculate the total distribution function hn (.g) of 
the whole channel. From Eq.(ll), hn{g) can be defined 
as: 



hnig) =jdai... d(7nf{(Ti) . . . f{(7n)5 



which, by using Eq.(14), reduces to: 



(16) 




(17) 



It is clear that hn{g) behaves as .g""" for g <C 1 since 
the integral appearing in the last equality of the above 
expression is well behaved in the g limit. In fact the 
g ^ limit of the Dirac (5- function appearing in Eq.(17) 
reduces to: 



lim 6 





(18) 



so that, finally: 



K{g) ^ (1 - a„)"'g /rf(7i ...do-„^(5(cri) rQai 

= n(l - a„)5~"" for g <^ 1. (19) 

The above equation is the main result of this paper, i. 
e., the distribution function of the occupied channel con- 
ductivities hn{g) is of the same form of Eq.(3). In this 
situation, for sufficiently large values of a„ the RRN con- 
ductivity can behave in a nonuniversal way with expo- 
nent t > 2.0. The condition for universality breakdown 
is given by Eq.(4) which for a three-dimensional network 
implies > ac — 0.107. From Eqs.(13,15) this condi- 
tion corresponds to: 



n <nc 



I - etc 
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(20) 



so that, for fixed values of and the value of 

the transport exponent t is governed solely by the num- 
ber of spheres that can be arranged within the occupied 
one-dimensional channels. The overall behavior of hn{g) 
is reported in Fig. 2 where we report a numerical calcula- 
tion of Eq.(16) (solid lines) together with the asymptotic 
behavior obtained in Eq.(19) (dotted lines). In this ex- 
ample we have set ^ = 2 nm, $ = 10 nm and L = 0.1 /itm 
corresponding to L/^ = 10, = 0.2 and Uc ~ 8.17. 
For n = 9 > ric the distribution function goes to zero 
as Eq.(19) with a„ ~ —0.244 while for n = 6 < ric 
hn{g) diverges for g — > with exponent a„ = 0.7. Since 
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FIG. 2. Distribution function hn{g) of the conductivity of 
the occupied channels for = 10, — 0.2 and different 
values of n. Solid lines are the result of a numerical calculation 
of Eq.(16) while the dotted lines are the asymptotic results of 
Eq.(19). 



Uc — 0.107, we expect that for n = 9 transport is univer- 
sal while for n = 6 the exponent t becomes larger than 
to ^ 2.0 as in Eq.(4). 

Before discussing our numerical results on the RRN 
conductivity, it is worth to point out that our model 
can be easily generalized to consider also situations in 
which the number of spheres accommodated in the one- 
dimensional channels is not fixc;cl. More spccdfically, if 
Vn' is the distribution function of the number n' of 
spheres, then the distribution function of the occupied 
channels is generalized to 



f^ia) = '^'Pn'hn'ig)- 



(21) 



As an instructive case let us consider a bimodal distribu- 
tion of the form: 



(22) 



where < q < 1. For q ~ or q = I wc recover the pre- 
vious case in which the occupied channels have the same 
number of spheres and whether transport is universal or 
not depends on the specific values nl (for g = 1) or n2 
(for (7=1). An interesting case is given by < g < 1 
and n2 < Uc < nl according to which there is a con- 
centration q of channels conductivities with distribution 
function with exponent a„i < Uc and 1 — q channels with 
cxn2 > etc- This case is depicted in Fig. 3 where nl = 9 
and n2 = 6 and, as in Fig. 2, L/$ = 10 and = 0.2. 
For all values q < 1 the 5 — > limit is governed by the 
diverging part of the total distribution function. In this 
case we expect that at the critical point the transport ex- 
ponent in not universal for any value q < 1. Note however 



FIG. 3. Channel conductivity distribution function h{g), 
Eq.(21), for the bimodal distribution of Eq.(22) and for dif- 
ferent values of q. L/^ = 10, ^/^ = 0.2, nl = 9 and n2 = 6. 



that for q sufficiently close to the unity, the asymptotic 
regime is reached for relatively small values of the c;on- 
ductivity. As we shall see in the next section, this has 
the effect of shrinking the region where criticality sets in 
with t = v — ane)- 



III. MONTE CARLO RESULTS ON THE CUBIC 
LATTICE 

In this section we discuss our Monte Carlo calculations 
for the conductivity a of the RRN model defined in the 
last section. In constructing the RRN we must first im- 
plement numerically the conductivity of the channels oc- 
cupied by a given number n of spheres. If x, (i = 1, . . . , n) 
is a set of random numbers equally distributed in the 
interval (0, 1) then it is easily found that the channel 
conductivity g having Eq.(17) as its corresponding dis- 
tribution function is: 



E 



(23) 



The RRN is then defined to have a fraction p of channels 
(in a cubic lattice) with g as given by Eq.(23) and a 
fraction 1—p with g = 0. The generalization to a bimodal 
distribution, Eqs.(21,22), is straightforward. 

To calculate numerically the transport exponent t we 
use the transfer-matrix method of Derrida and Vanni- 
menus applied to a simple cubic lattice of — 1 sites 
in the z direction. A'' sites along y and £ along the x 
direction. Periodic boundary conditions are used in the 
y-direction while to the top plane is applied a unitary 
voltage and the bottom plane is grounded to zero.^^'^^ 
For jC sufficiently large {£. N) this method permits 
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to calculate the conductivity per unit length of a cubic 
lattice. We calculate the conductivity for different 
linear sizes N at the percolation thresold Pc — 0.2488126 
for bond percolation on a cubic lattice, and then we 
extract by least-square fits the critical exponent t from 
the finite size scaling relation:^^'^^ 

CTjv = a7V-*/^(l + 6Af-'^), (24) 

where v ~ 0.88 is the correlation length exponent, a and 
h are constants and lo is the first correction to the scaling 
exponent t/v. In performing the calculations we have 
considered the following geometries: AT = 6 (£ = 5 x lO''), 
iV = 8 (£ = 2 X 10^), iV 10 (£ = 1 X 10^), TV = 12 
(£ = 8 X 106), iV = 14 (£ = 2 X 10^), and = 16 
(£ = 2 X 10^). 

In Fig. 4 we report the obtained values of the criti- 
cal exponent t for = 0.2 and for two different values 
of the ratio L/<^ between the length channel and sphere 
diameter. Each square corresponds to a particular num- 
ber n of inter-sphere tunneling junctions arranged in the 
channel (see the caption) which, from Eqs.(13,15), also 
gives the corresponding value of the tunneling exponent 
an reported in the abscissa. As a function of a;„, the criti- 
cal exponent t nicely follows Eq. (4) (solid curve) confirm- 
ing that universal {t ~ 2.0) or nonuniversal [t > 2.0) be- 
havior is obtained just by changing the number of spheres 
accommodated in the channels. In our least square fit- 
tings to Eq.(24) we have found that the minimum is 
obtained by setting b ^ and w ~ 1.0 for a„ < ttc and 
6 = for a„ > ttc- It is worth noticing that our Monte 
Carlo results on the cubic lattice agrees with Eq.(4) much 
better than the corresponding problem, Eqs.(2,3), on the 
two-dimensional square lattice. '^^''^^ 

We have applied the transfer-matrix method also to 
the bimodal distribution of Eq.(22) with L/^ = 10, 
= 0.2, nl = 9 and n2 = 6 and have found that, 
as expected, at Pc the critical exponent is nonuniversal 
already for q = 0.9. However what is interesting in the 
bimodal case is the behavior of the conductivity a away 
from the critical thresold. The high-structured shape 
of h{g) for < g < 1 reported in Fig. 3 in fact sug- 
gests that the p dependence of a could be affected by 
the competition between the two exponents q;„i < 
and a„2 > ac- To study this problem, the application 
of the transfer-matrix method for values of the occupied 
channel concentration p away from the critical thresold 
Pc is not efficient since the computational time of the al- 
gorithm increases as p is moved from pc-^^'^'^ Hence we 
have approached the problem by solving the RRN by the 
conjugate gradient method which is more efficient away 
from the critical point. ^'^ The resulting a is reported in 
Fig. 5 for a cubic lattice of 40 x 40 x 40 sites and periodic 
boundary conditions applied to the sides not connected 
with the external potential drop. We have considered 
the bimodal case defined by L/$ — 10, nl = 9, n2 = 6 
and different values of q. For q = 1 (filled squares in 
fig. 5) all the occupied channels have nl = 9 tunneling 
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FIG. 4. Critical exponent t as a function of the tunnel- 
ing exponent for = 0.2 and different values of L/$ 

and of the number n of inter-sphere tunneling junctions ac- 
commodated within the occupied channels of a cubic ran- 
dom-resistor network. From left to right: n = 9, 8, . . . , 5 for 
L/$ = 10 (filled squares) and n = 46, 45, 43, 41, 39, 33, 27 
for L/$ = 50 (open squares). The solid curve is the the- 
oretical result t = — 2.0 for a„ < Qc — 0.107 and 
1 = 1/ + 1/(1 - Q„) for an > otc [see Eq.(4)]. 

junctions and the conductivity is well approximated by 

Eq.(l) with critical exponent t = 1.8 ± 0.1. This value 
is slightly less than the universal result f = tg — 2.0 and 
this difference signals the limitation of extracting criti- 
cal exponents from the p dependence of a in finite size 
samples. However a follows the power-law Eq.(l) in the 
interval p — pc < 0.1 — 0.2. A nice power-law is found 
also for g = (filled diamonds) for which the occupied 
channels have n2 = 6 number of junctions. In this case 
however the exponent is t = 3.7 ± 0.2, i. e., slightly 
less than the nonuniversal value t ~ 4.0 obtained by the 
transfer- matrix method (see Fig. 4). 

In Fig. 5 we report also a calculated for intermediate 
values of q. For q = 0.9 and q = 0.8 the pc dependence 
of the conductivity can be reasonably fitted by a simple 
power-law only for p — Pc <~ 0.05 for which we have 
found t = 3.2± 0.3 and t = 4.1 ± 0.4, respectively. We 
interpret this shrinking of the critical region as being due 
to the large contribution of the fraction of channels with 
nl = 9 to the occupied channels distribution function 
h(g). However for q = 0.6 the critical region is already 
full restored and a follows a power-law with t = 3.8 ±0.5 
for p — Pc< 0.1. 

IV. DISCUSSION AND CONCLUSIONS 

As shown in the previous sections, the interesting char- 
acteristic of our segregated tunneling-percolation model 
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FIG. 5. Conductivity cr of a 40 x 40 x 40 cubic lattice with 
the bimodal distribution of Eq.(22) with nl = 9, n2 = 6, 
L/$ = 10 and = 0.2. Symbols are mean values of 10 dif- 
ferent runs with standard deviations given by the error bars. 
Dashed lines are fits to Eq.(l) with t = 1.8 ± 0.1 for g = 1 
and t = 3.7 ± 0.2 for g = 0. 



is the possibility of having universal or nonunivcrsal be- 
havior of transport within the same theoretical frame- 
work. As we have discussed, if the microscopic physical 
and geometric parameters (insulating CTibe size L, sphere 
diameter <D>, number of spheres and localization length ^) 
are such that the tunneling factor q,i is larger than the 
critical value Qc — 0.107 then the critical exponent t 
is nonuniversal and follows t = ^"-1-1/(1 — q;„), other- 
wise transport is universal and the critical exponent is 
t = to — 2.0. This universal/nonuniversal crossover is 
experimentally observed in thick-film resistors for which 
is reported to vary between t ~ 2.0 and t ~ 5.0 also for 
mixtures of chemically identical constituents. It can be 
argued that different fabrication procedures (for exam- 
ple firing temperature) affects the microstructure leading 
to different effective values of a„. Of course our model 
is oversimplified in the sense that interactions between 
the conducting and insulating phases are completely ne- 
glected. However it is remarkable that only two as- 
sumptions, quasi-one dimensionality of the conducting 
channels and Poisson distribution of the position of the 
spheres inside the channels, are sufficient to give rise to 
such rich phenomenology. 



The model discussed in this paper captures the es- 
sential physics, but eventually it can be further gener- 
alized to include more realistic features. For example, it 
is possible to account for different sizes of the conducting 
spheres in a straightforward manner, since also for this 
case the one-dimension nearest-neighbour distance dis- 
tribution function is provided by an analytical and exact 
expression. Also the tunneling expression Eq.(5) can 
be refined by including, for example, charging energies 
or distribution functions for the tunneling factor ^. 
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